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1 Introduction 

Ever since the initial applications of quantum mechanics to the dynamics of 
' electrons in solids, the analysis of Bloch electrons moving in a homogeneous 

, electric field has been of central importance. 

By employing semiclassical arguments, in 1928 Bloch demonstrated 
, that, a wave packet given by a superposition of single-band states peaked 

about some quasimomentum, ?ik, moves with a group velocity given by the 
gradient of the energy-band function with respect to the quasimomentum 
fH ' and that the rate of change of the quasimomentum is proportional to the 

Q I applied field F. This is often referred to as the "acceleration theorem": 

hk = eF . (1) 



o 

> 
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Thus, in the absence of interband tunneling and scattering processes, the 
quasimomentum of a Bloch electron in a homogeneous and static electric 
. field will be uniformly accelerated into the next Brillouin zone in a repeated- 



zone scheme (or equivalently undergoes an Umklapp process back in to the 
first zone). The corresponding motion of the Bloch electron through the pe- 
riodic energy-band structure, shown in Fig. |l|, is called "Bloch oscillation" ; 
It is characterized by an oscillation period tb = h/eFd, where d denotes 
the lattice periodicity in the field direction. 

There are two mechanisms impeding a fully periodic motion: interband 
tunneling and scattering processes. Interband tunneling is an intricate prob- 
lem and still at the center of a continuing debate. Early calculations of the 
tunneling probability into other bands in which the electric field is rep- 
resented by a time-independent scalar potential were made by Zener |^ 
using a Wentzel-Kramers-Brillouin generalization of Bloch functions, by 
Houston using accelerated Bloch states (Houston states), and subse- 
quently by Kane Q and Argyres |^ who employed the crystal-momentum 
representation. Their calculations lead to the conclusion that the tunnel- 
ing rate per Bloch period is much less than unity for electric fields up to 
lO^V/cm for typical band parameters corresponding to elemental or com- 
pound semiconductors. 

Despite the apparent agreement among these calculations, the validity 
of employing the crystal-momentum representation or Houston functions 
to describe electrons moving in a non-periodic (crystal plus external field) 
potential has been disputed. The starting point of the controversy was the 
original paper by Wannier He pointed out that, due to the transla- 
tional symmetry of the crystal potential, if (j){r) is an eigenfunction of the 
scalar-potential Hamiltonian (corresponding to the perfect crystal plus the 
external field) with eigenvalue e, then any (j){r+nd) is also an eigenfunction 
with eigenvalue e -I- riAe, where Ae = eFd is the so-called Wannier-Stark 
splitting (d being the primitive lattice vector along the field direction) . He 
concluded that the translational symmetry of the crystal gives rise to a 
discrete energy spectrum, the so-called Wannier-Stark ladder. The states 
corresponding to these equidistantly spaced levels arc localized states, as 
schematically shown in Fig. || for the case of a semiconductor superlattice. 
The degree of this Wannier-Stark localization depends on the strength of 
the applied field. 

The existence of such energy quantization was disputed by Zak , who 
pointed out that for the case of an infinite crystal the scalar potential — F-r 
is not bounded, which implies a continuous energy spectrum. Thus, the 
main point of the controversy was related to the existence (or absence) of 
Wannier-Stark ladders. More precisely, the point was to decide if interband 
tunneling (neglected in the original calculation by Wannier ) is so strong 
to destroy the Wannier-Stark energy quantization (and the corresponding 
Bloch oscillations) or not. We will give a brief historical account of this 
long-standing controversy in the following section. 

It is only during the last decade that this controversy came to an end. 
From a theoretical point of view, most of the formal problems related to 
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the non-periodic nature of the scalar potential (superimposed to the pe- 
riodic crystal potential) were finally removed by using a vector-potential 
representation of the applied field ^. Within such vector-potential pic- 
ture, upper boundaries for the interband tunneling probability have been 
established at a rigorous level, which show that an electron may execute a 
number of Bloch oscillations before tunneling out of the band ^ |ll| , 
in qualitatively good agreement with the earlier predictions of Zener and 
Kane[||. 

The second mechanism impeding a fully periodic motion is scattering by 
phonons, impurities, etc. (see Fig. This results in lifetimes shorter than 
the Bloch period tb for all reasonable values of the electric field, so that 
Bloch oscillations should not be observable in conventional solids. 

In superlattices, however, the situation is much more favourable because 
of the smaller Bloch period tb resulting from the small width of the mini- 
Brillouin zone in the field direction [p^ . 

Indeed, the existence of Wannier-Stark ladders as well as Bloch oscilla- 
tions in superlattices has been confirmed by a number of recent experiments 
[p^ . The photoluminescence and photocurrent measurements of the biased 
GaAs/GaAlAs superlattices performed by Mendez and coworkers 14 1, to- 
gether with the electroluminescence experiments by Voisin and cowork- 
ers fist , provided the earliest experimental evidence of the field- induced 
Wannier-Stark ladders in superlattices. A few years later, Feldmann and 
coworkers [|l^ were able to measure Bloch oscillations in the time domain 
through a four-wave- mixing experiment originally suggested by von Plessen 
and Thomas ||l^ . A detailed analysis of the Bloch oscillations in the four- 
wave-mixing signal (which reflects the interband dynamics) has been also 
performed by Leo and coworkers |l^ . 

In addition to the above interband-polarization analysis, Bloch oscilla- 
tions have been also detected by monitoring the intraband polarization 
which, in turn, is reflected by anisotropic changes in the refractive index 
[ p^ . Measurements based on transmittive electrooptic sampling (TEOS) 
have been performed by Dekorsy and coworkers |2^, |2l[. Finally Bloch 
oscillations were recently measured through a direct detection of the Ter- 
aHertz (THz) radiation in semiconductor superlattices |2|, [2^ . 

The aim of this chapter is to present a general approach to the study of 
the ultrafast carrier dynamics in semiconductor superlattices. Our theoreti- 
cal description, based on the density matrix formalism discussed in chapter 
6, is presented in Sect. It allows to derive a set of kinetic equations which 
accounts for interband tunneling as well as scattering processes and it is 
valid in any quantum-mechanical representation. 

In Sect, m the freedom of choice of the basis states in our kinetic for- 
mulation will be used to introduce the two typical pictures commonly 
used for the description of semiconductor superlattices, namely, the Bloch- 
oscillation and Wannier-Stark pictures. In particular, we will see that they 
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correspond to the two equivalent vector- and scalar-potential representa- 
tions of the applied field. This will implicitly state the total equivalence 
of the Bloch-oscillation and Wannier-Stark representations, which in turn 
shows that the so-called "semiclassical Bloch picture" is on the contrary a 
rigorous quantum-mechanical result. 

Finally, in Sect. ^ we will review and discuss some simulated experiments. 

2 Hysterical background 

In this section we give a brief historical account concerning the contro- 
versy on the existence of Wannier-Stark ladders mentioned above. Some of 
the main criticisms to the pioneering works on Bloch oscillations may be 
summarized as follows. 

(i) The eigenvalues of the time-independent Schrodinger equation are not 
quantized but they form a continuous spectrum. 

(ii) Since the Hamiltonian within the scalar-potential representation is not 
periodic, it is not clear whether one can employ the periodic Bloch states 
or Houston functions: a superposition of Bloch functions will automati- 
cally yield a periodic function while the solution of the time-dependent 
Schrodinger equation is, in general, not periodic. 

(iii) The position operator (entering the scalar potential) is hill-defined with- 
in the crystal-momentum representation. 

In a series of papers, Wannier Q and coworkers |2^, have argued 
that in the presence of a homogeneous electric field, one can modify the 
Bloch states in such a way that there is no interband coupling and an 
electron in a crystal will move within one band with its k changing in time 
according to the acceleration theorem in Eq. (^. Furthermore, if k(i — 0) 
is in the direction of a reciprocal-lattice vector, the periodic motion in 
k-space gives rise to an energy quantization with Ae — eFd {d = ^ 
being the lattice constant along the field direction) , the so-called Wannier- 
Stark ladders. "The basis for this idea is that energy bands arise from the 
translational symmetry of the crystalline field and this symmetry is not 
removed physically by the presence of the applied field" (24|, i.e. the field is 
still periodic with the lattice period. 

These arguments have been refuted by Zak 0, who shows that, al- 
though it immediately follows from the one-dimensional time-independent 
Schrodinger equation for an infinite crystal with lattice constant d, that 
if e is an eigenvalue, e + nAe is also an eigenvalue, the spectrum of e is 
continuous with — oo < e < +oo, so the Wannier-Stark ladders do not 
exist. 

Wannier has argued that Zak's criticisms of his proof are not valid, 
but concedes that the Stark ladders may be metastable resonant states 
limited by interband tunneling, as for the case of the hydrogen atom in 
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the presence of a static field. However, Wannier's arguments were imme- 
diately rejected by Zak |2^, who claims that Wannier's original equation 
was incorrect. 

A few years later, Rabinovitch and Zak have extended Zak's ear- 
lier arguments to the question of Bloch oscillations. They argue that since 
neglecting the interband coupling terms in the time-independent Schrodin- 
ger equation leads incorrectly to energy quantization, then the interband 
terms cannot be neglected in the lowest approximation because they are 
the same order as the terms retained. By applying the same reasoning to 
the time-dependent equation, they conclude (without offering a proof) that 
neglecting the interband terms as a first approximation, as done by Hous- 
ton is incorrect for times equal to or longer than the Bloch-oscillation 
period tb ■ From their conclusions it follows that the typical diagrams which 
are commonly used to portray trajectories of k(i) superimposed upon the 
energy-band structure (see Fig. ^ are incorrect and misleading. 

Nevertheless, shortly before these latter arguments appeared, experimen- 
tal results were obtained by Koss and Lambert , which were interpreted 
as supporting the existence of Wannier-Stark levels. They found that the 
observed low-temperature optical absorption of GaAs in a strong electric 
field {F = lO^V/cm) closely followed the theoretical predictions of Call- 
away which, in turn, were based on employing Kane's wavefunctions 
and Wannier-Stark quantized energy levels. 

3 Theoretical analysis 

In this section we will try to review and discuss in a systematic way the 
basic ideas used in the theoretical analysis of semiconductor superlattices. 
As already pointed out in Sect. the phenomena under investigation, i.e. 
Bloch oscillations and Wannier-Stark localization, are peculiar of any lattice 
structure. Therefore, even if most of the results discussed in this chapter 
refer to semiconductor superlattices, the general formulation presented in 
this section applies to any crystalline structure. 

3.1 Physical system 

In order to study the optical and transport properties of semiconductor 
superlattices, let us consider a gas of carriers in a crystal under the action of 
an applied electromagnetic field. The carriers will experience their mutual 
interaction as well as the interaction with the phonon modes of the crystal. 
Such physical system can be described by the following Hamiltonian: 

H = He + Hp + Hcc + Hep -t- Hpp ■ (2) 

The first term describes the noninteracting-carrier system in the pres- 
ence of the external electromagnetic field while the second one refers to 
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the free-phonon system. The last three terms describe many-body con- 
tributions: they refer, respectively, to carrier-carrier, carrier-phonon, and 
phonon-phonon interactions. 

In order to discuss their explicit form, let us introduce the usual second- 
quantization field operators 'J'^^i-) and *(r). They describe, respectively, 
the creation and the annihilation of a carrier in r. In terms of the above 
field operators the carrier Hamiltonian He can be written as 



2rno 



. (3) 



Here, V^(v) denotes the periodic potential due to the perfect crystal while 
A(r,t) and iy9(r, i) denote, respectively, the vector and scalar potentials 
corresponding to the external electromagnetic field. Since we are interested 
in the electrooptical properties as well as in the ultrafast dynamics of pho- 
toexcited carriers, the electromagnetic field acting on the crystal — and the 
corresponding electromagnetic potentials — will be the sum of two different 
contributions: the high-frequency laser field responsible for the ultrafast 
optical excitation and the additional electromagnetic field acting on the 
photoexcited carriers on a longer time-scale. More specifically, by denoting 
with the labels 1 and 2 these two contributions, we can write 

A(r,i) = Ai(r,t)+A2(r,t) , ^{v ,t) ^ ^^{y ,t) ^ ^^{v ,t) (4) 

and recalling that 

E(r,i) = -i|-A(r,f)-Vr(^(r,t) , B(r, i) - Vr x A(r, i) (5) 

c at 

we have 

E(r,i) =Ei(r,i)+E2(r,i) , B(r, t) = Bi(r, i) + B2(r, i) . (6) 

Equation , which gives the electromagnetic fields in terms of the corre- 
sponding vector and scalar potentials, reflects the well known gauge free- 
dom: there is an infinite number of possible combinations of A and Lp 
which give rise to the same electromagnetic field {E, B}. We will use such 
freedom of choice for the laser field (term 1): we assume a homogeneous 
(space-independent) laser field Ei(i) fully described by the scalar potential 

(/9i(r,<) = -Ei(t)-r . (7) 

This assumption, which corresponds to the well known dipole approxima- 
tion, is well justified as long as the space-scale of interest is small compared 
to the light wavelength. The explicit form of the laser field considered in 
this chapter is 
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where Eo{t) is the ampUtude of the hght field and ul denotes its central 
frequency. 

With this particular choice of the electromagnetic potentials describing 
the laser field, the Hamiltonian in Eq. (||) can be rewritten as 

H, = + H,, , (9) 



where 

> 2 



(-zWr-fA2(r,<))' 



2mc 



+ eip2{r,t) + V'{r) 



*(r) (10) 



describes the carrier system in the crystal under the action of the electro- 
magnetic field 2 only, while 

H,, = ey'dr*t(r)^i(r,i)*(r) (11) 

describes the carrier-light (cl) interaction due to the laser photoexcitation. 

In analogy with the carrier system, by denoting with 5^ ^ and the 
creation and annihilation operators for a phonon of mode A and wavevector 
q, the free-phonon Hamiltonian takes the form 

Hp^E^'^q^^qA^A , (12) 

qA 

where ujqx is the dispersion relation for the phonon mode A. 

Let us now discuss the explicit form of the many-body contributions. 
The carrier-carrier interaction is described by the two-body Hamiltonian 

Hec = i Jdr y"dr'*t(r)*t(i.')i/^(r-r')*(r')*(r) , (13) 

where Vcc denotes the Coulomb potential. 

Let us now introduce the carrier-phonon interaction Hamiltonian 

- I dr*t(r)l/,p(r)*(r) , (14) 

where 

- E [sqA^Ae^'' '- + rcxbl.e-^-^-^] (15) 

qA 

is the electrostatic phonon potential induced by the lattice vibrations. Here, 
the explicit form of the coupling function depends on the particular 
phonon mode A (acoustical, optical, etc.) as well as on the coupling mech- 
anism considered (deformation potential, polar coupling, etc.). 

Let us finally discuss the phonon-phonon contribution Upp. The free- 
phonon Hamiltonian Hp introduced in Eq. |l^, which describes a system 
of noninteracting phonons, by definition accounts only for the harmonic 
part of the lattice potential. However, non-harmonic contributions of the 
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interatomic potential can play an important role in determining the lattice 
dynamics in highly excited systems , since they are responsible for the 
decay of optical phonons into phonons of lower frequency. In our second- 
quantization picture, these non-harmonic contributions can be described in 
terms of a phonon-phonon interaction which induces, in general, transitions 
between free-phonon states. Here, we will not discuss the explicit form 
of the phonon-phonon Hamiltonian Hpp responsible for such a decay. We 
will simply assume that such phonon-phonon interaction is so efficient to 
maintain the phonon system in thermal equilibrium. This corresponds to 
neglect hot-phonon effects . 

It is well known that the coordinate representation used so far is not the 
most appropriate one in describing the electron dynamics within a peri- 
odic crystal. On the contrary, it is in general more convenient to employ 
the representation given by the eigenstates of the noninteracting-carrier 
Hamiltonian — or a part of it — since it automatically accounts for some 
of the symmetries of the system. For the moment we will simply consider 
an orthonormal basis set {(/)„(r)} without specifying which part of the 
Hamiltonian is diagonal in such representation. This will allow us to write 
down equations valid in any quantum-mechanical representation. Since the 
noninteracting-carrier Hamiltonian is, in general, a function of time, also 
the basis functions may be time-dependent. Here, the label n denotes, 
in general, a set of discrete and/or continuous quantmn numbers. In the 
absence of electromagnetic field, the above wavefunctions will correspond 
to the well known Bloch states of the crystal and the index n will reduce to 
the wavevector k plus the band index v. In the presence of a homogeneous 
magnetic field the eigenfunctions (/>„ may instead correspond to Landau 
states. Finally, for the case of a constant and homogeneous electric field, 
there exist two equivalent representations: the accelerated Bloch states and 
the Wannier-Stark picture. Such equivalence results to be of crucial im- 
portance in understanding the relationship between Bloch oscillations and 
Wannier-Stark localization and, for this reason, it will be discussed in more 
detail in Sect. |4[ 

Let us now reconsider the system Hamiltonian introduced so far in terms 
of such (j)n representation. As a starting point, we may expand the second- 
quantization field operators in terms of the new wavefunctions: 

*(r) = ^</.„(r)a„ , *t(r) = ^CWa], • (16) 

n n 

The above expansion defines the new set of second-quantization operators 
a\ and a„; They describe, respectively, the creation and annihilation of a 
carrier in state n. 

For the case of a semiconductor crystal (which will be the only one 
considered in this chapter), the energy spectrum of the noninteracting- 
carrier Hamiltonian in Eq. (|l^) — or a part of it — is always characterized 
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by two well-separated energy regions: the valence and the conduction band. 
Also in the presence of an applied electromagnetic field, the periodic lattice 
potential V' gives rise to a large energy gap. Therefore, we deal with two 
energetically well-separated regions, which suggests the introduction of the 
so-called electron-hole picture. This corresponds to a separation of the set 
of states {(pn} into conduction states {tp^} and valence states {(p'j}- Thus, 
also the creation (annihilation) operators ajj (a„) introduced in Eq. (^6|) 
will be divided into creation (annihilation) electron and hole operators: c| 
(Cj) and {dj). In terms of the new electron-hole picture, the expansion 
in Eq. ( p^ ) is given by: 

*(r) = ^0,^(r)c,+^0f(r)d] 



*t(r) = ^C(r)c,U^0^'(r)d,. 



(17) 



If we now insert the above expansion into Eq. ([l0|), the noninteracting- 
carrier Hamiltonian takes the form 



where 



jj' 



(18) 



2mo 



eip2 + V-eo 



<j>\l\v) (19) 



are just the matrix elements of the Hamiltonian in the 0-representation. 
The ± sign refers, respectively, to electrons and holes while Eq denotes the 
conduction-band edge. Here, we neglect any valence-to-conduction band 
coupling due to the external electromagnetic field and vice versa. This is 
well fulfilled for the systems and field-regimes we are going to discuss in 
this chapter. As already pointed out, the above Hamiltonian may be time- 
dependent. We will discuss this aspect in the following section, where we 
will derive our set of kinetic equations. 

Let us now write in terms of our electron-hole representation the carrier- 
light interaction Hamiltonian introduced in Eq. (llh: 



H 



cl — 



E 



-{t)d^c. 



(20) 



The above expression has been obtained within the well known rotating- 
wave approximation by neglecting intraband transitions, absent for the case 
of optical excitations. Here 
between states cpl and (j)^. 



fjtl^ denotes the optical dipole matrix element 
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Similarly, the carrier-carrier Hamiltonian ([L3h can be rewritten as 



, rt 

2 



- E (21) 



*1«2J1J2 



where 



^??2;3;4 = / / dr>r, (r)0r2 (r')^^^(r - r')^Z3 {^')^U (r) (22) 



are the Coulomb matrix elements within our (/)-representation. The first 
two terms describe the repulsive electron-electron and hole-hole interaction 
while the last one describes the attractive electron-hole interaction. Here, 
we neglect terms that do not conserve the number of electron-hole pairs, i.e. 
impact-ionization and Auger-recombination processes |Q, as well as the 
interband exchange interaction. This monopole-monopole approximation is 
justified as long as the exciton binding energy (which in semiconductors is 
less than 20meV) is small compared to the energy gap (which is more than 
leV). 

Finally, let us rewrite the carrier-phonon interaction Hamiltonian intro- 
duced in Eq. (|l^): 



Hep — ^ [sii'.qACi^qACi' + 5'ii*,qAcI'&qA'^i 
ii' ,qA 

E [5',',qA4V^.'+ff^*,qArf,UlArf,] (23) 



with 



at^^ = ~9^>.\dv<i>t\r)e''-^'lH!\r). (24) 

In Eq. (^) we can clearly recognize four different contributions correspond- 
ing to electron and hole phonon absorption and emission. 



3.2 Kinetic description 

Our kinetic description of the ultrafast carrier dynamics in semiconductor 
superlattices, presented in this section, is based on the density-matrix for- 
malism. Since this approach has been already reviewed and discussed in 
chapter 6, here we will simply recall in our notation the kinetic equations 
relevant for the analysis of carrier dynamics in superlattices, generaliz- 
ing the approach of chapter 6 to the case of a time-dependent quantum- 
mechanical representation. 
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The set of kinetic variables is the same considered in chapter 6. Given 
our electron-hole representation {(pf}, {4>'j}: we will consider the intraband 
electron and hole single-particle density matrices 



fu'^{4c^') , flr={d]dr) (25) 

as well as the corresponding interband density matrix 

Pji = {djC,) . (26) 

Here, the diagonal elements ff^ and /j^ correspond to the electron and 
hole distribution functions of the Boltzmann theory while the non-diagonal 
terms describe intraband polarizations. On the contrary, the interband 
density-matrix elements p^^ describe interband (or optical) polarizations. 

In order to derive the set of kinetic equations, i.e. the equations of motion 
for the above kinetic variables, the standard procedure starts by deriving 
the equations of motion for the electron and hole operators introduced in 
Eq. (0): 

q = j dr0r(r)*(r) , = j dr^f (r)*t(r) . (27) 

By applying the Heisenberg equation of motion for the field operator 'S', 



I.e. 



it is easy to obtain the following equations of motion: 



r^. = AK'H]+iE^«'^/ (29) 



with 



Z:/^^^hJdr(^l^f'*{r)yl/\r). (30) 

As for the case of Eq. (p^, here we neglect again valence-to-conduction 
band coupling and vice versa. Compared to the more conventional Heisen- 
berg equations of motion, the above equations contain an extra-term, the 
last one. It accounts for the possible time dependence of our (/)-representa- 
tion which will induce transitions between different states according to the 
matrix elements Zn, . 

By combining the above equations of motion with the definitions of 
the kinetic variables in Eqs. (E5[]29), our set of kinetic equations can be 
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schematically written as: 





H 




H 



d_ 




d 




d 



(31) 



H 



They exhibit the same structure of the equations of motion (^9|) for the 
electron and hole creation and annihilation operators: a first term induced 
by the system Hamiltonian H (which does not account for the time varia- 
tion of the basis states) and a second one induced by the time dependence 
of the basis functions ((>. 

Let us start discussing this second term, whose explicit form is given by: 



Such terms were not considered in chapter 6, where a time-independent 
representation has been used. As we will see in Sect. they will play 
a central role for the description of Bloch oscillations within the vector- 
potential representation. 

Let us now come to the first term. This, in turn, is the sum of different 
contributions, corresponding to the various parts of the Hamiltonian. In 
particular, the total Hamiltonian can be regarded as the sum of two terms, 
a single-particle contribution plus a many-body one: 



H = Hsp + Hmb — (H° + Hcj + Hp) + (Hcc + Hep + Hpp) . (33) 



The explicit form of the time evolution due to the single-particle Hamil- 
tonian Hsp (non-interacting carriers plus carrier-light interaction plus free 
phonons) is given by: 






(32) 
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This is a closed set of equations, which is a consequence of the single- 
particle nature of H^p. In addition, we stress that the structure of the 
two contributions entering Eq. ( ^l] ) is very similar: one can include the 
contribution (|3^) into Eq. ( ^4| ) by replacing e with e + Z. 

Let us finally discuss the contributions due to the many-body part of the 
Hamiltonian: carrier-carrier and carrier-phonon interactions (the phonon- 
phonon one is not explicitly considered here). As discussed in chapter 6, 
for both interaction mechanisms one can derive a hierarchy of equations in- 
volving higher-order density matrices and, in order to close such equations 
with respect to our set of kinetic variables, approximations are needed. 
The lowest-order contributions to our equations of motion are given by 
first-order terms in the many-body Hamiltonian: Hartree-Fock level. Since 
we will neglect coherent-phonon states, the only Hartree-Fock contribu- 
tions will come from carrier-carrier interaction. They simply result in a 
renormalization 

of the single-particle energy matrices e*^/'' as well as in a renormalization 
Af/n,i=-E^^5u...P.... (36) 

i2j2 

of the external field U. (The explicit form of the renormalization terms 
considered in this chapter accounts for the Fock contributions only, i.e. 
no Hartree terms. The general structure of Hartree-Fock contributions, 
relevant for the case of a strongly non-homogeneous system, is discussed in 
chapter 6.) 

We stress that the Hartree-Fock approximation, which consists in fac- 
torizing average values of four-point operators into products of two density 
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matrices, is independent from the quantum-mechanical picture. This is a 
general property: any mean-field approximation gives the same result in 
different representations. The reason is that the mean- field operation com- 
mutes with any unitary transformation connecting different basis states. It 
is then clear that the above kinetic equations are valid in any quantum- 
mechanical representation. 

All the contributions to the system dynamics discussed so far describe a 
fully coherent dynamics, i.e. no scattering processes. In order to treat inco- 
herent phenomena, e.g. energy relaxation and dephasing, one has to go one 
step further in the perturbation expansion taking into account also second- 
order contributions (in the perturbation Hamiltonian Umb)- The deriva- 
tion of these higher-order contributions, discussed in chapter 6, will not 
be repeated here. Again, as for the first-order contributions (Hartree-Fock 
terms), in order to obtain a closed set of equations (with respect to our set 
of kinetic variables ( p5| - |2^ )) additional approximations are needed, namely 
the mean-field and the Markov approximation. As for the Hartree-Fock 
case, the mean-field approximation allows to write the various higher-order 
density matrices as products of single-particle ones. The Markov approx- 
imation allows to eliminate the additional higher-order kinetic variables, 
e.g. phonon-assisted density matrices, providing a closed set of equations 
still local in time, i.e. no memory effects |3^, ^ This last ap- 
proximation is not performed in the quantum-kinetic theory discussed in 
chapter 6 where, in addition to our single-particle variables, one considers 
two-particle and phonon-assisted density matrices ||^, ^ . 

While the mean- field approximation is representation-independent, this 
is unfortunately not the case for the Markov limit. This clearly implies 
that the validity of the Markov approximation is strictly related to the 
quantum-mechanical representation considered. We will come back to this 
point in the following section where the two different pictures used for the 
study of the carrier dynamics in superlattices are discussed. 

The above kinetic description, based on intra- and interband density 
matrices, allows us to evaluate any single-particle quantity. In particular, 
for the analysis of semiconductor superlattices two physical quantities play 
a central role: the intra- and interband total (or macroscopic) polarizations: 

P'^/'^it) = J2 Kl'^fin^^) . P^'^^Y. : (37) 

W ij 

where M'^^'^ and ^"^^ denote, respectively, the intra- and interband dipole 
matrix elements in our (/)-representation. The time derivative of the intra- 
band polarization P^/'* describes the radiation field induced by the Bloch- 
oscillation dynamics (which for a superlattice structure is in the TeraHertz 
range) while the Fourier transform of the interband (or optical) polarization 
peh TpYOvides the optical-absorption spectrum. 
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4 Two equivalent pictures 



In this section we will apply the theoretical approach presented so far to the 
case of a semiconductor superlattice in the presence of an uniform (space- 
independent) electric field. The non-interacting carriers within the super- 
lattice crystal will then be described by the Hamiltonian H° in Eq. (p^, 
where now the electrodynamic potentials A2 and kpi (in the following sim- 
ply denoted with A and 93) correspond to a homogeneous electric field 
E2(r,0=F(i). 

As pointed out in Sect. |3.l| , the natural quantum- mechanical represen- 
tation is given by the eigenstates of this Hamiltonian: 

'(-»Wr-fA(r,t))' 



2mo 



e(^(r,i) + y'(r) 



</)„(r) = e„</)„(r) . (38) 



However, due to the gauge freedom discussed in Sect. 3T, there is an infinite 
number of possible combinations of A and — and therefore of possible 
Hamiltonians — which describe the same homogeneous electric field F(i). In 
particular, one can identify two independent choices: the vector-potential 
gauge ^ 

A(r,i) = -c\ F{t')dt' , Lp{r,t) = (39) 

and the scalar-potential gauge 

A(r,i)=0, ^(r,i)--F(t)-r (40) 

(previously employed for the description of the laser photoexcitation in 
Eq. (|)). 

As we will see, the two independent choices correspond, respectively, to 
the well known Bloch-oscillation and Wannier-Stark pictures. They simply 
reflect two equivalent quantum-mechanical representations and, therefore, 
any physical phenomenon can be described in both pictures. 



The Bloch-oscillation picture 



The vector-potential approach presented in this section, originally proposed 
by Kittel 1^, is discussed in Ref. 

Within the vector-potential gauge (^9|), the above eigenvalue equation 
reduces to: 



:A(t)y 



2mo 



V\r) 



0n(r) = en(?!'ri(r) 



(41) 



In this gauge the vector potential is space- independent but, even for the 
case of a static field (F(t) = Fo), it is always time-dependent. Therefore, 
the above Hamiltonian (together with its eigenvalues and eigenfunctions) 
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will be time dependent as well. However, for any time t we can consider 
its "instantaneous" eigenstates 0„(r, i). They can be easily evaluated by 
means of the following transformation: 

0„(r,i)=C(r,t)e^^("^-*) (42) 

with 

X(r,t)-A(0T . (43) 

By applying this transformation to the eigenvalue problem in Eq. (^), we 
obtain |9|: 



2mo 



(r) = e„C(r) , (44) 



i.e. the wavefunctions f/)" (r) are just the Bloch states ?!'^iy(r) of our semi- 
conductor, and the energy spectrum e„ coincides with the carrier band 
structure eki^ Therefore, from Eq. ( p^ ) the desired eigenfunctions result to 
be of the form: 

0k.(r,i) = 0L(r,i)e^^(*'-'- . (45) 
Apart from a phase-factor, they coincide with the conventional Bloch states 
of the crystal. The reason can be understood as follows: Also in the presence 
of the electric field F, the Hamiltonian in Eq. ( ^ ) is still invariant under 
a lattice translation corresponding to the crystal potential Thus, the 
crystal momentum fik is still a "good" quantum number and the band 
dispersion remains the same: e = eki^ Therefore, at each time t our time- 
dependent eigenstates seem to coincide (a part from the phase-factor) with 
the Bloch states of the crystal and, at a first glance, it is not clear which 
is the role played by the applied field. 

In order to answer this question, let us consider again the general form 
of our eigenstates in Eq. (^). At the initial time to the vector potential 
A is equal to zero and, therefore, the two basis sets coincide: '/>k^j,(r, to) = 
(/)^j,(r). (Here, ko and k denote the carrier wavevectors at time to and 
t, respectively. They are, in principle, independent quantities, since they 
correspond to two different eigenvalue problems.) In other words, the Bloch 
states can be also regarded as the states (f)^^^^ at the initial time to and 
vice versa. This allows us to rewrite Eq. (Ea) as 



0k.(r,t) = 0k<,.(r,to)eSAW-r (45) 

or equivalently 

KM = C.We-^^'*^-'- . (47) 
Moreover, in view of the translational symmetry of the crystal, at each time 
t the Bloch states (t>° should obey the Bloch theorem ||] : 

(/)£,(r + a) = 0°,(r)e^'^-, (48) 

where a denotes any periodicity vector of the crystal. If we now apply the 
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Bloch theorem to both sides of Eq. (|47|), i.e. at time t and to, we finaUy 
obtain: 

k = ko - ^A(t) . (49) 
nc 

This result is quite important: the symmetry properties of the Hamiltonian 
require a precise relationship between the (formally independent) wavevec- 
tors ko and k. More specifically, the carrier wavevector results to be a 
function of time (k = k(t)), i.e. the instantaneous sets of basis states 
{0k(t)i'(r, i)} (corresponding to different times t) are mutually connected 
through a continuous time evolution of the wavevector k. We can then an- 
swer the previous question saying that within the vector-potential gauge 
the application of a homogeneous field F{t) induces a simple "drift" in k- 
space of the crystal Bloch states, which are therefore not "distorted" by 
the presence of the field. 

From a physical point of view, the above equation describes the contin- 
uous time evolution of a carrier in k-space induced by the applied electric 
field. In particular, taking into account the explicit form of the vector po- 
tential A given in Eq. (^9|), we have 

k(t) = ko + J / F{t')dt' , (50) 

from which the acceleration theorem in Eq. (|^) is recovered: 

m = • (51) 

As pointed out in the introduction, this is usually regarded as a "semiclas- 
sical" result, i.e. as obtained by applying to a Bloch electron the laws of 
classical mechanics. However, the above analysis shows that this is a rigor- 
ous quantum- mechanical result: the quantum evolution of a carrier within 
a given band under the action of a homogeneous electric field is rigorously 
described by the acceleration theorem. 

We want to stress once again that within the vector-potential gauge 
discussed so far the acceleration theorem is just a result of the symmetry 
properties of the crystal and the time-dependent eigenstates in Eq. ( |45[ ) 
describe the quantum analogue of the "semiclassical motion" of a carrier 
in k-space (see Fig. 0). 

Let us now consider the case of a static field, i.e. F{t) — Fo, applied 
parallel to a symmetry axis of the crystal. In this case, the corresponding 
vector potential entering Eq. (^ ) is a linear function of time, which induces 
a uniform drift of the carriers in k-space along the field direction: 

k(i) =ko+k(t-to) = ^(t-io) . (52) 



Since the carrier energy — given by the eigenvalue in Eq. mW) — coincides 
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with the crystal band structure eki^, its time evolution is: 

e,{t) = e,{k{t)) = eJko+k{t-to)) . (53) 



Due to the periodic nature of the band structure as a function of k, i.e. 
ey(k) = ey(k + G) (G being a reciprocal-lattice vector), the carrier will 
execute a periodic motion in time with a period 

where d is the lattice periodicity (in real space) along the field direction. 
This coincides with the Bloch period previously introduced (see Sect. |l|). 
It corresponds to the time needed for the electron to travel from any point 
k to the energetically equivalent point fc + ^ . 

These periodic oscillations of the carrier over the crystal band structure 
(see Fig. |l]) are known as Bloch oscillations. As pointed out in Sect. 0, 
they were first introduced by Bloch ^ on the basis of semiclassical argu- 
ments. However, as for the case of the acceleration theorem discussed above, 
this is a rigorous quantum-mechanical result of the vector-potential pic- 
ture discussed so far. Such a clear physical interpretation of the quantum- 
mechanical theory in terms of a semiclassical picture is hard to obtain 
within the scalar-potential gauge presented in the following section. 

Both the acceleration theorem ( |5l|) and the Bloch-oscillation dynamics 
previously discussed are induced by the noninteracting-carrier Hamiltonian 
in Eq. ( ^l|) through its time-dependent eigenstates (j)k(t)iy Therefore, the 
Bloch-oscillation dynamics considered so far does not account for many- 
body effects (carrier-carrier and carrier-phonon interactions) as well as for 
the effects induced by the time variation of our basis states. For a more 
"realistic" description of the carrier dynamics within our vector-potential 
picture we are then forced to employ the general kinetic theory presented 



in Sect. 3.2 



As discussed in chapter 6, for the case of a homogeneous semiconductor 
crystal the only relevant terms of the single-particle density matrix in our 
hp representation are those diagonal in k. This property, which is due again 
to the translational symmetry of the Hamiltonian, reduces the set of kinetic 
variables in Eqs. (psfpq) to the intraband density-matrix elements 



/k,Qa' — y'^ka'^ka'/ ' /-k,/3/3' ^ \'^-k/3'^-k/3' / ^^^^ 

plus the interband density-matrix elements 

Pk,l3a = (c'-k/jCka) • (56) 



Here, the standard electron-hole picture introduced in Sect. 3.1 has been 
applied to our set of time-dependent eigenstates 0k(t)i': the band index v 
(which refers to both conduction and valence states) is replaced by two 
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separate band indices a and P for electrons and holes, respectively, while, 
due to the charge-conjugation symmetry, the hole states are still labeled in 
terms of the corresponding valence-electron states, i.e. k''/3 = — k*^/?. 

Let us now discuss the explicit form of the kinetic equations (|^) in our 
vector-potential picture: 



re re 



H 



d 

di 



re 

J k,Qia2 



rh _ ^ rh 



k,/3i/32 



H 



d 
It 



a rh 

"ir/-k,/3i/32 



d 

+ ^Pk,ftai 



H 



(57) 



Here, the first term is induced by the system Hamiltonian H — H^p + Umb 
(which does not account for the time variation of the basis states) while 
the second one is induced by the time evolution of the basis functions 4>. 

The contributions to the carrier dynamics due to the single-particle 
Hamiltonian Hgp are given in Eq. (Q). They consist in a "free rotation" 
plus a term due to the interaction with the external laser field. For the case 
of an ultrafast laser excitation, after the initial carrier photogeneration the 
only non-vanishing contributions in Eq. (Q) are such free-rotation terms. 
If we now consider that in our vector-potential representation the energy 
matrix e in Eq. (n9) is diagonal, 



e/h 
^k;,k'/' 



ek('''^kk'<5;i' 



(58) 



the single-particle contributions after the initial photoexcitation reduce to: 



±re 



[^ka2 ^kaj /k,aiQ2 



sp 



dt 



rh 



132 



TT [^-k/32 <^-k/3i] /-k,ft/32 



d 



1 

ih 



-k/3i 



^kaj Pk,/3iai 



(59) 



As we can see, the above equations describe a set of independent many- 
level systems, i.e. one for each k value. In addition, there is no coupling 
between different kinetic variables. If, in particular, we assume a diagonal 
initial condition 



re/h 



e/h c- 



(60) 
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(which is well fulfilled for the case of a laser photoexcitation of standard 
bulk semiconductors as well as superlattice structures), the kinetic equa- 
tions ( p^ for the non-zero density-matrix elements, i.e. for the diagonal 
ones, reduce to: 

(01) 

where the ± sign refers, respectively, to electrons and holes. If we now 
remember that in our vector-potential representation the wavevector k is 
itself a function of time (see Eq. (^9|)), the above kinetic equations can be 
rewritten as: 

|/±k' ± k • = , + k • VkPk. - . (62) 

For both distribution functions / and interband polarizations p we obtain 
a simple drift equation whose general solution is of the form: 

y(k(t),i) = j/(ko,io) , (63) 

i.e. the function at time t is obtained through a rigid shift Ak = k — ko 
of the hmction at the initial time to- Such drift in k-space, induced by the 
external field F, is again described by the acceleration theorem in Eq. (^l|). 
Therefore, as expected, the carrier dynamics described by the above ki- 
netic equations (which accounts for the single-particle Hamiltonian only) 
is fully equivalent to the Bloch-oscillation dynamics discussed above. How- 
ever, such relatively simple picture of the carrier motion does not account 
for the time-dependence of our basis states as well as for many-body effects, 
e.g. carrier-carrier and carrier-phonon interactions. 

The contributions to the carrier dynamics induced by the time variation 
of the basis states are given in Eq. (^2|). In our vector-potential represen- 
tation ( ^ ) the explicit form of the matrix elements Z^^,^ introduced in 
Eq. (|^ is 

^kfk'. = <«"5ku' (64) 

with 

Z^J, = ±e {5u' - 1) 1 dr<^^/'" [Y{t) ■ r] 0^/," . (65) 

They result to be strictly related to the matrix elements of the scalar po- 
tential in Eq. (^o|), as we will discuss in the following section. The k/ k'l' 
transitions induced by the time variation of the basis states are always di- 
agonal in k; this reflects the momentum conservation in the carrier-field in- 
teraction, i.e. since the momentum q corresponding to a space-independent 
field is equal to zero, the initial and final carrier wavevectors coincide. More- 
over, there are no intraband (l = I') transitions, which confirms that the 
action of the field within a given band is fully described by the drift terms 
in Eq. (|6|). 

If we now rewrite Eq. (^) in our vector-potential representation taking 
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into account the explicit form of the above matrix elements, we finally 
obtain: 



dt 



re 

J k,ai(: 



d rh 



''k,/32/34^/3i/33 ^-k,/33/3i'^/32/34] /-k.fe/Ji 



dt 



1 

ih 



/33/34 

E 

Q2/32 



.■^-k,/3i,32'^"ia2 



(66) 



From the above kinetic equations we clearly see that the time variation 
of our basis states (j> induces "vertical" (i.e. k = k') transitions between 
different bands. Such interband coupling is the well known Zener tunneling 
This effect is usually described as an interband transition induced by 
the scalar potential — F • r, which is also evident from the explicit form 
of the Zener matrix elements in Eq. (|65|). However, within our vector- 
potential picture the Zener tunneling originates from the time variation 
of our accelerated Bloch states in Eq. (|4^). 

As discussed in Ref. |^ , for the case of bulk semiconductors this interband 
coupling results to be very limited even for the case of high applied fields. 
Therefore, the Bloch-oscillation scenario of the semiclassical theory (see 
Fig. |l|) is practically unmodified by Zener tunneling. For the case of inter- 
est, i.e. that of semiconductor superlattices, interminiband Zener tunneling 
is expected to play a significant role in the high-field regime. However, for 
relatively low fields (up to lO^V/cm) the effect is again negligible and the 
Bloch-oscillation regime is fully recovered. In this case, the time-scale of 
Zener-tunneling processes is much longer than the Bloch-oscillation period 
tb- Therefore, the effect due to the time variation of our basis states is 
negligible, i.e. the time variation can be regarded as an "adiabatic trans- 
formation" . 

Let us finally consider the role played by many-body effects, i.e. carrier- 
carrier and carrier-phonon interactions. As discussed in Sect. 3.2 as well as 
in chapter 6, these many-body effects can be divided into coherent and in- 
coherent contributions. With coherent contributions we refer to first-order 
terms in the many-body Hamiltonian Umb- Since we neglect coherent- 
phonon states, the only non-zero contributions originate from carrier-carrier 
interaction. The explicit form of these Hartree-Fock terms in our vector- 
potential representation is obtained from Eqs. ( |35[]3^ ) by replacing the 
generic labels i and j with ka and — k/3, respectively. 

From a physical point of view, these coherent contributions give rise to 
excitonic and band-renormalization effects, which result in a modification of 
the single-particle energy spectrum. In our case, these excitonic effects may 
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lead to modifications of the Bloch-oscillation dynamics, e.g. small variations 
of the Bloch period. As we will see, within the Wannier-Stark picture such 
excitonic effects manifest themselves in a modification of the single-particle 
Wannier-Stark energy levels, which are not equidistantly spaced anymore 

Let us now come to the role played by incoherent contributions, i.e. 
second-order contributions in the many-body Hamiltonian. As pointed out 
in Sect. 3.2, these terms are usually treated within the usual Markov ap- 
proximation. The resulting contributions describe, in general, second-order 
transitions connecting all possible kinetic variables, i.e. all possible density- 
matrix elements. The transitions connecting diagonal density-matrix el- 
ements, i.e. distribution functions, can be easily described in terms of 
stochastic scattering processes, i.e. due to the scattering with a partner car- 
rier or with a phonon, the electron may undergo a transition from an initial 
state ki/ to a final state kV. On the contrary, for the transitions involving 
non-diagonal terms the second-order coupling is not positive-definite, i.e. 
it is not a rate, and the intuitive scattering picture cannot be employed. 

It is not in the spirit of this chapter to derive and discuss the explicit form 
of the second-order carrier-carrier and carrier-phonon contributions. From 
a physical point of view, both carrier-carrier and carrier-phonon scattering 
processes give rise to energy relaxation and dephasing. It is well known that, 
due to scattering events, a photogenerated carrier distribution will relax 
both energy and momentum and, in addition, it will lose its internal degree 
of coherence, which corresponds to a decay of the interband polarizations. 
This stochastic dynamics may strongly influence the deterministic Bloch- 
oscillation regime discussed so far [^0[ ^ 

As we will see from some simulated experiments reported in Sect. ||, the 
role played by carrier-carrier and carrier-phonon scattering strongly de- 
pends on the physical conditions considered, e.g. carrier density, excitation 
energy, and lattice temperature. When the scattering rate corresponding to 
the dominant interaction mechanism is much larger than the Bloch oscilla- 
tion frequency = the Bloch oscillations are not suppressed, i.e. the 
carriers perform on average several Bloch oscillations between two scatter- 
ing events. On the contrary, if the scattering rate is larger than the Bloch 
frequency, the carrier cannot execute a full oscillation without scattering. 
In this case, the Bloch oscillations are totally suppressed and we deal with 
a diffusive-transport regime. 

As discussed in the introductory part of this chapter, in bulk semicon- 
ductors also for very high fields the Bloch-oscillation period is larger than 
the typical scattering times. On the contrary, in semiconductor superlat- 
tices the Bloch period is at least one order of magnitude smaller than in 
bulk systems and, therefore, comparable or even smaller than the typical 
scattering times. 

This allows us to answer the controversial question: "do Bloch oscillations 
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really exist?". The analysis presented in this section shows that, in the 
absence of scattering events, Bloch oscillations exist and, contrary to the 



early papers by Rabinovitch and Zak 1 29 , they are not significantly affected 



by Zener tunneling, both for bulk and super lattices. 

On the contrary, due to scattering events, Bloch oscillations are fully 
suppressed in bulk semiconductors but they still survive in superlattices, 
as confirmed by several experiments ||l^, |l^, |2^, ^ . 

In the following section we will discuss the so-called Wannier-Stark pic- 
ture. Contrary to the vector-potential approach discussed so far, this will 
correspond to a scalar-potential representation. In particular, we will study 
the link between the two pictures showing that phenomena which are pe- 
culiar of one picture can be equally described within the second one. 



4-2 The Wannier-Stark picture 

Within the scalar-potential gauge (^), the eigenvalue equation (^) reduces 
to: 

fi^V^ 1 

" y'(r)-eFT 0(r) =e</.(r) . (67) 



2too 

In this gauge the scalar potential is space dependent but, for the case of a 
static field considered in this section, it is always time independent. 

As originally pointed out by Wannier M, due to the translational sym- 
metry of the crystal potential (r) — V^{r + d), d being the primitive 
lattice vector along the field direction, if 0(r) is an eigenfunction of the 
Hamiltonian in Eq. ( |67| ) with eigenvalue e, then any (/)(r 4- nd) is also an 
eigenfunction with eigenvalue e -t- nAe, where Ae = eFd is the so-called 
Wannier-Stark splitting. However, as pointed out by Zak |^], for the case of 
an infinite crystal the scalar potential — F • r is not bounded, which implies 
a continuous energy spectrum. This was the starting point of the long- 
standing controversy on the existence of Wannier-Stark ladders discussed 
in Sect. ||. 

Today, after three decades, we know that the problem is by itself ill- 
defined and eigenstates can be defined only asymptotically pO| . In particu- 
lar, there exist ladders of metastable Wannier-Stark states weakly coupled 
(through Zener tunneling) to a continuous energy spectrum. 

We will now study the explicit form of these Wannier-Stark states. As 
a starting point, let us write the generic eigenstate (f) as a superposition 
of conventional Bloch states cj)'^^, i.e. let us move to the so-called crystal- 
momentum representation (CMR): 

1 f^^ 

</'W / c?ksk+G^0k+G,.(r) , (68) 

ki/ Gi/ 

where G is a generic reciprocal-lattice vector while D, denotes the volume 
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of the first Brillouin zone. Due to the periodicity of the Bloch states in k- 
space ((/"k+Gi/ ~ '/'ki^)' allowed to limit the above expansion to the 

first Brillouin zone, which implies to impose the same periodicity on the 
coefficients Ski^: Sk+Ci/ — Skz^- By inserting the above expansion (limited to 
the first Brillouin zone) in Eq. (|67|), our eigenvalue problem can be rewritten 
as: 



y- 



c^k'iJki/.kV'Sk'i.' = eski. , (69) 



where 

^kiy.k'i/' — eki/^ki/,k'i^' + ^kL,k'i/' C^^) 

are the matrix elements of the Hamiltonian H° in the CMR. The first term 
corresponds to the perfect (field-free) crystal while the second one describes 
the scalar-potential term: 

ffk.,k'.' = -e 1 drcj^Z [F . r] . (71) 

Following the approach discussed in Ref. [Q, the above scalar-potential 
matrix elements can be divided into intraband {v — v') and interband 
(y 7^ v') terms; The intraband terms can always be written as a drift 
operator @ || Q 

1 /"^ii d 
- / dk[Hl,^^,, = -iei^4,fe' ^ (72) 



(fc|| and fcx being, respectively, the components of the wavevector k parallel 
and perpendicular to the field) while the non-diagonal terms coincide with 
the Zener-tunneling matrix elements Z in Eq. (p3) . As already pointed out 



in Sect. 4.1, the Zener tunneling, which in the vector-potential represen- 
tation is induced by the time variation of the basis states, corresponds to 
interband transitions induced by the scalar potential in Eq. (|4C|). 

For moderate values of the applied field, Zener tunneling to other bands 
can be neglected and, by inserting Eqs. ( [70|j7^ ) into Eq. (|6^), our eigenvalue 
equation reduces to 

3 i 

-^Ski/ = ^ (ek,y - e) Ski/ , (73) 

whose solution is given by: 

Sk. = 5fc^fc^e II . (74) 

As expected, the coefficients are diagonal with respect to k±, i.e. the linear 
combination in Eq. ( |68| ) will only involve Bloch states with the same per- 
pendicular component k±. This reflects the translational symmetry of the 
Hamiltonian with respect to the plane perpendicular to the field. Moreover, 
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the periodicity condition Sk+c v — Ski^ (apphed along the field direction) 
requires that 

—^l (efc||fc^i. - e) = 27rn (75) 

which, in turn, tells us that the only allowed energy values are 

e = e^-L""^ = efc-LO-- + n eFd (76) 

with 

What we obtain is the Wannier-Stark ladder mentioned above, whose cen- 
tral {n = 0) value is given by the average of the band along the field 
direction (for a given k±). 

Let us now discuss the corresponding wavefunctions. By inserting in the 



expansion (68) the explicit form of the coefficients Skz^ given in Eq. (|7J) we 
have: 

^'--(r) A |^^'&||e^^^''^.^--^^^"^^'^'=V^„,^.(r) . (78) 

They are the so-called Wannier-Stark states shown in Fig. |[ For each band 
u (and for any given fcj_), we have a set of energetically equidistant states 
^k±nv ^ obtained from the central (n — 0) state (f)^^^'^ through the spatial 
translation r v+nd discussed above. As schematically depicted in Fig. ||, 
each state is localized around one of the atomic cells of the superlattice 
and the degree of localization depends on the strength of the applied field 
F. More precisely, when the Wannier-Stark energy eFd is much smaller 
than the superlattice miniband width, the wavefunction (j)^-^™ are weakly 
localized, they extend over several elementary cells. On the contrary, when 
eFd is comparable or larger than the miniband width, the localization 
increases and the function results to be significantly different from zero only 
in one cell. Since the miniband width of the holes is smaller than that of 
the electrons, the hole states exhibit a stronger Wannier-Stark localization 
(see Fig. |). 

Due to the neglect of interband Zener tunneling, each Wannier-Stark 
state (f)^^™ in Eq. ( [T^ ) is obtained as a linear combination of Bloch states 
belonging to the same band v. Moreover, we see that all Bloch states have 
the same weight in the expansion, i.e. the coefficients are just phase-factors. 

The coefficients s introduced so far can be regarded as the matrix ele- 
ments of a unitary transformation connecting the Bloch to the Wannier- 
Stark representation. It is then clear that the inverse transformation allows 
us to write any Bloch state as a linear combination of Wannier-Stark states: 
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The Wannier-Stark states in Eq. ( [Tq ) can now be used as basis states 
for our kinetic description. The kinetic variables in the Wannier-Stark rep- 
resentation wiU be formaUy the same as for the vector-potential picture 



discussed in Sect. 4.1. They are defined according to Eqs. ( |55| , p6| ) where 
the three-dimensional wavevector k is replaced by its perpendicular com- 
ponent k±^ while the band index a//3 is replaced by the same band index 
plus the Wannier-Stark ladder index n. 

Provided the above label substitution, the time evolution of the new ki- 
netic variables is again described by Eq. (|57|). However, within our Wannier- 
Stark representation the basis states are time-independent and, therefore, 
the 0-contributions in Eq. (|5^) are equal to zero. Moreover, contrary to the 
vector-potential case, the energy matrix e in Eq. ( [l9|) is not diagonal. The 
diagonal terms are now given by the Wannier-Stark ladders in Eq. ( [76|) 
while the non-diagonal terms are once again Zener-tunneling matrix ele- 
ments between different Wannier-Stark states. 



As already pointed out in Sect. 3.2, the coherent contributions entering 



our kinetic equations are independent from the quantum-mechanical repre- 
sentation considered. This tells us that the Zener-tunneling contributions to 
the equations of motion in the scalar- and vector-potential representations 
should coincide. In fact, in spite of their different physical interpretations 
(in the vector-potential picture they are induced by the time variation of 
the basis states while in the scalar-potential one they are due to inter- 
band transitions induced by the field Hamiltonian), their formal structure 
is exactly the same, as can be seen by comparing Eqs. (^) and (|4|). 

On the contrary, incoherent contributions, i.e. scattering terms derived 
within Markov approximation, are expected to be representation-depen- 
dent. In particular, within the vector-potential picture the scattering terms 
are usually derived by neglecting the time dependence of the basis states 
0k(t)i/i which corresponds to neglecting the action of the field, i.e. the time 
variation of the carrier wavevector k(t), during the collision. Within the 
Wannier-Stark picture, the basis states are time-independent and the stan- 
dard Markov limit automatically accounts for the so-called intracollisional 
field effect |3^, However, for moderate values of the applied field 
incoherent contributions evaluated in the scalar- and vector-potential rep- 
resentations coincide. 

Before concluding this section, let us try to better understand the physi- 
cal link between the accelerated Bloch states of the vector-potential picture 
and the Wannier-Stark states of the scalar-potential representation. Both 
basis sets have been introduced as eigenstatcs of two equivalent Hamil- 
tonians, corresponding to the two different electromagnetic gauges (see 
Eqs. (^) and (^7|)). However, it is well known |^ that the solutions of 
the corresponding time-dependent Schrodinger equations coincide (a part 
from a phase- factor which is physically irrelevant). More specifically, let us 
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consider the generic time-dependent Schrodinger equation 

ihj^i^{v,t)^H{t)^{v,t) (80) 

together with the corresponding eigenvalue problem 

H{t)MT^,t)^€x{t)Mv,t) . (81) 

The generic solution ip at time t is given by a linear combination of the 
eigenstates (j) according to the time evolution induced by the Hamiltonian 
H: 

V'(r,t) =^5Ae"^^*^^'*'^'*'(/.A(r,t) , (82) 

A 

which, for the case of the vector-potential Hamiltonian (^TJ) and the corre- 
sponding eigenstates in Eq. ( |45| ) reduces to: 

^(r,i) =e'I^^W--^5k„.e-'^/'o'''<"'^'*'C(,),(r) , (83) 

where ko = k(io) denotes the carrier wavevector at the initial time to- As 
discussed in Ref . , it is always possible to consider the gauge transforma- 
tion connecting the above vector-potential picture to the scalar-potential 
one. The generator of this gauge transformation is the function —A{t) ■ r, 
which tells us that, going from the vector- to the scalar-potential picture, 
the first phase-factor in Eq. ( |83| ) cancels exactly with the corresponding 
phase-factor of the gauge transformation. Thus, the time-dependent wave 
function (|3|) in the scalar-potential gauge reduces to: 

i,{v, t)^Yl ^ko.e"'^ '^'^''"'UmAr) . (84) 

This is a linear combination of the so-called Houston states 1^ originally 
introduced as time-dependent solutions of the scalar-potential Schrodinger 
equation. 

On the other hand, within the Wannier-Stark representation ([78|), the 
linear combination in Eq. ( ^^ reduces to: 

V;(r,t) = ^fc^n^e-^<^'-'""(*-*°V'^"''(i-) . (85) 

It is then clear that for a given initial condition ■!/)(r, to) the two last linear 
combinations must give at any time t the same wavefunction ^. Let us 
consider as initial condition a single Wannier-Stark state (j)^-^^'^ . According 
to Eq. (^5|), at time t the function ip differs from that at time to only by 
a phase-factor, which implies that this will be a stationary state, i.e. the 
wavefunction will remain localized around a given cell and \ip\'^ will not 
change in time. Therefore, it should be possible to choose the coefficients 5 
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entering Eq. ( p4[ ) in such a way that the corresponding expansion in terms 
of Houston states will provide the same stationary Wannier-Stark state. 

In order to determine the explicit form of the coefhcients corresponding to 
a stationary state, let us rewrite Eq. (|8^), replacing at each time t the sum 
over ko with an equivalent sum over the instantaneous k = ko + k(i — to) 
(see Eq. (H)): 

=E^k-k(*-*„).e^"^-^^'""'''''^L(r) . (86) 

The stationary-state condition corresponds to impose that each individual 
term entering the above expansion will evolve in time according to the 
constant Wannier-Stark energy e'^-L"'^^ 

5,_,(,_,„),e-^^>^<-'"'^*'(xe-i^^-"(*-*»), (87) 
which implies that for each time t 

The above time integral over t' can be translated into a corresponding 
integral over k' — k\\{t'): 

Since the first phase- factor on the right-hand side is time-independent, we 
finally obtain: 

5k„.oce^^C(^-'"-^''""'')'"=' . (90) 

As expected, the coefficients 5 of a stationary state coincide with the co- 
efficients s in Eq. ([z^). Thus, the linear combinations of accelerated Bloch 
states corresponding to stationary states are just the Wannier-Stark states 
introduced in Eq. ( |7q ) as eigenstates of the scalar-potential Hamiltonian. 

From the above analysis, we see that within a time-dependent approach 
the scalar- and vector-potential pictures are totally equivalent. According 
to the initial condition (i.e. depending on the coefficients S), we may have a 
Bloch-oscillation as well as a Wannier-Stark scenario, or any intermediate 
regime. If we consider, for example, the case of a laser excitation whose 
energy spectrum is concentrated around a well defined frequency, this will 
generate a distribution of photoexcited carriers strongly peaked about a 
particular k. Each carrier will then be described by a single accelerated 
Bloch state (pki, and will execute Bloch oscillations. Thus, the overall mo- 
tion of this packet in k-space will resemble the periodic motion of a single 
electron "prepared" in a Bloch state (see Fig. ^). As we will discuss in the 
following section, such Bloch-oscillation dynamics can be monitored via 
four-wave-mixing experiments or THz-signal measurements. 
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On the contrary, if we perform a linear-absorption measurement using 
as a light source a laser with uniform spectral distribution, we generate a 
uniform distribution of photoexcited carriers in k-space with their proper 
phase coherence (described by the corresponding interband polarizations). 
Such distribution is shifted in k-space but, being almost constant, there 
is no macroscopic effect. This resembles the situation corresponding to a 
single electron prepared in a Wannier-Stark state, which is a uniform su- 
perposition of Bloch states. This is confirmed by optical-absorption inves- 
tigations which clearly show the Wannier-Stark energy quantization (see 
Fig. |). 



5 Some simulated experiments 

In this section, we will review recent simulated experiments of the ultra- 
fast carrier dynamics in semiconductor superlattices ^ ^ ^ . 
They are based on a generalized Monte Carlo solution |5C|, |52|, |33[ of the 
set of kinetic equations (so-called semiconductor Bloch equations) derived 



in Sect. 3.2. This generalized Monte Carlo approach, successfully applied 
for the interpretation of ultrafast coherent phenomena in bulk semiconduc- 
tors |5^ , is based on a combined solution of our kinetic equations 
the coherent contributions are evaluated by means of a direct numer- 
ical integration while the incoherent ones are "sampled" by means of a 
conventional Monte Carlo simulation in the three-dimensional k-space. 
This generalized Monte Carlo method has been recently applied to semi- 



conductor superlattices. As described in Refs. |41, 46|, the simulation sche- 



me is based on the Bloch-state representation of the vector-potential picture 



introduced in Sec. 4.1 



The following superlattice model has been employed: The energy disper- 
sion and the corresponding wavefunctions along the growth direction (fc||) 
are computed within the well known Kronig-Penney model p^ , while for 
the in-plane direction (fc_L) an effective- mass model has been used. Starting 
from these three-dimensional wavefunctions (j)^^, the various carrier-carrier 
as well as carrier-phonon matrix elements are numerically computed (see 
Eqs. ( p2|) and (^J)). They are, in general, functions of the various mini- 
band indices and depend separately on k^^ and k±, thus fully reflecting the 
anisotropic nature of the superlattice structure. 

Only coupling to GaAs bulk phonons has been considered. This, of 
course, is a simplifying approximation which neglects any superlattice ef- 
fect on the phonon dispersion, such as confinement of optical modes in the 



wells and in the barriers, and the presence of interface modes [58 . How- 
ever, while these modifications have important consequences for phonon 
spectroscopies (like Raman scattering), they are far less decisive for trans- 
port phenomena. Indeed, by now it is well known M, B9| that the total 
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scattering rates are sufficiently well reproduced if the phonon spectrum is 
assumed to be bulk-like. 

We will start discussing the scattering-induced damping of Bloch oscilla- 
tions. In particular, we will show that in the low-density limit this damping 
is mainly determined by optical-phonon scattering ^] while at high 
densities the main mechanism responsible for the suppression of Bloch os- 
cillations is found to be carrier-carrier scattering . 

This Bloch-oscillation analysis in the time domain is also confirmed by 
its counterpart in the frequency domain. As pointed out in Sect. ^, the pres- 
ence of Bloch oscillations, due to a negligible scattering dynamics, should 
correspond to Wannier-Stark energy quantization. This is confirmed by 
the simulated optical-absorption spectra, which clearly show t he p resence 
of the field-induced Wannier-Stark ladders introduced in Sect. [4.2| [^ . 



5.1 Bloch-oscillation analysis 

All the simulated experiments presented in this section refer to the su- 
perlattice structure considered in Ref. |Q: 111 A GaAs wells and 17 A 
Alo.3Gao.7As barriers. For such a structure there has been experimental 
evidence for a THz-emission from Bloch oscillations |^ . 

In the first set of simulated experiments an initial distribution of pho- 
toexcited carriers (electron-hole pairs) is generated by a 100 fs Gaussian 
laser pulse in resonance with the first-miniband exciton {Ulol ~ 1540 meV). 
The strength of the applied electric field is assumed to be 4 kV/cni, which 
corresponds to a Bloch period tb = h/eFd of about 800 fs. 

In the low-density limit (corresponding to a weak laser excitation), in- 
coherent scattering processes do not alter the Bloch-oscillation dynamics. 
This is due to the following reasons: In agreement with recent experimental 



[ p3[ , 40 and theoretical |4l|, ^ investigations, for supcrlattices char- 
acterized by a miniband width smaller than the LO-phonon energy — as 
for the structure considered here — and for laser excitations close to the 
band gap, at low temperature carrier-phonon scattering is not permitted. 
Moreover, in this low-density regime carrier-carrier scattering plays no role: 
Due to the quasi-elastic nature of Coulomb collisions, in the low-density 
limit the majority of the scattering processes is characterized by a very 
small momentum transfer; As a consequence, the momentum relaxation 
along the growth direction is negligible. As a result, on this picosecond 
time-scale the carrier system exhibits a coherent Bloch-oscillation dynam- 
ics, i.e. a negligible scattering-induced dephasing. This can be clearly seen 
from the time evolution of the carrier distribution as a function of k\\ (i.e. 
averaged over fcj^) shown in Fig. ||. During the laser photoexcitation [t = 0) 
the carriers are generated around k\\ = 0, where the transitions are close to 
resonance with the laser excitation. According to the acceleration theorem 
(P, the electrons are then shifted in k-space. When the carriers reach the 



30 



border of the first Brillouin zone they are Bragg reflected. After about 800 
fs, corresponding to the Bloch period r^, the carriers have completed one 
oscillation in k-space. As expected, the carriers execute Bloch oscillations 
without loosing the synchronism of their motion by scattering. This is again 
shown in Fig. || where we have plotted: (b) the mean kinetic energy, (c) 
the current, and (d) its time derivative which is proportional to the emitted 
far field, i.e. the THz-radiation. (It can be shown that, by neglecting Zener 
tunneling, the intraband polarization P'^/'' in Eq. ( p^ is proportional to 
the current.) All these three quantities exhibit oscillations characterized by 
the same Bloch period Tg. Due to the finite width of the carrier distri- 
bution in k-space (see Fig. §(a)), the amplitude of the oscillations of the 
kinetic energy is somewhat smaller than the miniband width. Since for this 
excitation condition the scattering-induced dephasing is negligible, the os- 
cillations of the current are symmetric around zero, which implies that the 
time average of the current is equal to zero, i.e. no dissipation. 

As already pointed out, this ideal Bloch-oscillation regime is typical of 
a laser excitation close to gap in the low-density limit. Let us now dis- 
cuss, still at low densities, the case of a laser photoexcitation high in the 
band. Figure ||(a) shows the THz-signal as obtained from a set of sim- 
ulated experiments corresponding to different laser excitations |46| . The 
different traces correspond to the emitted THz-signal for increasing exci- 
tation energies. We clearly notice the presence of Bloch oscillations in all 
cases. However, the oscillation amplitude and decay (effective damping) is 
excitation-dependent . 

For the case of a laser excitation resonant with the first- miniband exciton 
considered above (see Fig. ||), we have a strong THz-signal. The amplitude 
of the signal decreases when the excitation energy is increased. Additionally, 
there are also some small changes in the phase of the oscillations, which 
are induced by the electron-LO phonon scattering. 

When the laser energy comes into resonance with the transitions between 
the second electron and hole minibands {hujL ~ 1625 meV), the amplitude 
of the THz-signal increases again. The corresponding THz-transients show 
an initial part, which is strongly damped and some oscillations for longer 
times that are much less damped. For a better understanding of these 
results, we show in Fig. ^(b) the individual THz-signals, originating from 
the two electron and two heavy-hole minibands for the excitation with 
hu! = 1640 meV. The Bloch oscillations performed by the electrons within 
the second miniband are strongly damped due to intra- and intcrminiband 
LO-phonon scattering processes Since the width of this second 

miniband (45 meV) is somewhat larger than the LO-phonon energy, also 
intraminiband scattering is possible, whenever the electrons are accelerated 
into the high-energy region of the miniband. The THz-signal originating 
from electrons within the first miniband shows an oscillatory behavior with 
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a small amplitude and a phase which is determined by the time the electrons 
need to relax down to the bottom of the band. 

At the same time, the holes in both minibands exhibit undamped Bloch 
oscillations, since the minibands are so close in energy that for these ex- 
citation conditions no LO-phonon emission can occur. The analysis shows 
that at early times the THz-signal is mainly determined by the electrons 
within the second miniband. At later times the observed signal is due to 
heavy holes and electrons within the first miniband. 

The above theoretical analysis closely resembles experimental observa- 
tions obtained for a superlattice structure very similar to the one modelled 
here ||2^ . In these experiments, evidence for THz-emission from Bloch oscil- 
lations has been reported. For some excitation conditions these oscillations 
are associated with resonant excitation of the second miniband. The general 
behavior of the magnitude of the signals, the oscillations and the damping 
are close to the results shown in Fig. IJ. 

Finally, in order to study the density dependence of the Bloch-oscillation 
damping, let us go back to the case of laser excitations close to gap. Fig- 
ure ||(a) shows the total (electrons plus holes) THz-radiation as a function 
of time for three different carrier densities. With increasing carrier den- 
sity, carrier-carrier scattering becomes more and more important: Due to 
Coulomb screening, the momentum transfer in a carrier-carrier scattering 
increases (its typical value being comparable with the screening wavevec- 
tor). This can be seen in Fig. ||(a), where for increasing carrier densities 
we realize an increasing damping of the THz-signal. However, also for the 
highest carrier density considered here we deal with a damping time of the 
order of 700 fs, which is much larger than the typical dephasing time, i.e. 
the decay time of the interband polarization, associated with carrier-carrier 
scattering. The dephasing time is typically investigated by means of four- 
wave-mixing (FWM) measurements and such multi-pulse experiments can 
be simulated as well 1 54 , 5^ . From a theoretical point of view, a qualitative 
estimate of the dephasing time is given by the decay time of the "incoher- 
ently summed" polarization (ISP) |5l[. Figure ^(b) shows such ISP as a 
function of time for the same three carrier densities of Fig. 1(a). As ex- 
pected, the decay times are always much smaller than the corresponding 
damping times of the THz-signals (note the different time-scale in Figs, ^(a) 
and (b)). This difference can be understood as follows: The fast decay times 
of Fig. |^(b) reflect the interband dephasing, i.e. the sum of the electron and 
hole scattering rates. In particular, for the Coulomb interaction this means 
the sum of electron-electron, electron-hole, and hole-hole scattering. This 
last contribution is known to dominate and determines the dephasing time- 
scale. On the other hand, the total THz-radiation in Fig. ^(a) is the sum of 
the electron and hole contributions. However, due to the small value of the 
hole miniband width compared with the electron one, the electron contri- 
bution will dominate. This is clearly shown in Fig. ^, where the electron (a) 
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and hole (b) contributions to the THz-radiation are plotted as a function of 
time (note the different vertical scale). This means that the THz damping 
in Fig. ^(a) mainly reflects the damping of the electron contribution (see 
Fig. ||(a)). This decay, in turn, reflects the intraband dephasing of electrons 
which is due to electron-electron and electron-hole scattering only, i.e. no 
hole-hole contributions. This clearly explains the different decay times of 
Fig.|(a) andFig.|(b). 

From the above analysis we can conclude that the decay time of the 
THz-radiation due to carrier-carrier scattering differs considerably from 
the corresponding dephasing times obtained from a FWM experiment: The 
first one is a measurement of the intraband dephasing while the second one 
reflects the interband dephasing. 



5.2 Optical-absorption analysis 

Let us now discuss the frequency-domain counterpart of the Bloch-oscilla- 
tion picture considered so far. Similar to what happens in the time domain, 
for sufficiently high electric fields, i.e. when the Bloch period tb — h/eFd 
becomes smaller than the dephasing time, the optical spectra of the su- 
perlattice are expected to exhibit the frequency-domain counterpart of the 
Bloch oscillations, i.e. the Wannier-Stark energy quantization discussed in 



Sect. 4.2. In absence of Coulomb interaction, the Wannier-Stark ladder ab- 
sorption increases as a function of the photon energy in a step-like fashion. 
These steps are equidistantly spaced. This spacing, named Wannier-Stark 
splitting, is proportional to the applied electric field (see Eq. ([T^)). 

The simulated linear-absorption spectra corresponding to a superlattice 
structure with 95 A GaAs wells and 15 A Alo.3Gao.7As barriers are shown 
in Fig. 1^ |48j. As we can see, the Coulomb interaction gives rise to ex- 
citonic peaks in the absorption spectra and introduces couplings between 
these Wannier-Stark states. Such exciton peaks, which arc not equidistantly 
spaced any more, are often referred to as excitonic Wannier-Stark ladders 
[ p9[ of the superlattice. 

Since for the superlattice structure considered in this simulated experi- 
ment Q the combined miniband width is larger than the typical two- 
and three-dimensional exciton binding energies, it is possible to investigate 
the quasi-three dimensional absorption behavior of the delocalized mini- 
band states as well as localization effects induced by the electric field. 

For the free-field case, the electron and hole states are completely delo- 
calized in our three-dimensional k-spacc. The perturbation induced by the 
application of a low field (here « 5 kV/cm), couples the states along the 
field direction and in the spectra the Franz-Keldysh effect, well known from 
bulk materials, appears: one clearly notices oscillations which increase in 
amplitude with the field and shift with F^^^ from the n = and n = 1 
levels toward the center of the combined miniband. 
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For increasing field the potential drop over the distance of a few quantum 
wells eventually exceeds the miniband width and the electronic states be- 
come more and more localized. Despite the field-induced energy difference 
neFd, the superlattice potential is equal for quantum wells separated by 
nd. Therefore, the spectra decouple into a series of peaks corresponding to 
the excitonic ground states of the individual electron-hole Wannier-Stark 
levels. Each Wannier-Stark transition contributes to the absorption with a 
pronounced 1 — s exciton peak, plus higher bound exciton and continuum 
states. The oscillator strength of a transition n is proportional to the over- 
lap between electron and hole wavefunctions centered at quantum wells n' 
and n + n', respectively. The analysis shows that this oscillator strength is 
almost exclusively determined by the amplitude of the electron wavefunc- 
tion in quantum well n' + n since for fields in the Wannier-Stark regime the 
hole wavefunctions are almost completely localized over one quantum well 
due to their high effective mass (see Fig. ||). Thus, the oscillator strengths 
of transitions to higher \n\ become smaller with increasing \n\ and field. 

At high fields (here >w 8 kV/cm) the separation between the peaks is 
almost equal to neFd. For example, the peak of the n — transition which 
is shifted by the Wannier-Stark exciton binding energy with respect to the 
center of the combined miniband, demonstrates that the increasing local- 
ization also increases the exciton binding energy. This increased excitonic 
binding refiects the gradual transition from a three- to a two-dimensional 
behavior. 

For intermediate fields there is an interplay between the Wannier-Stark 
and the Franz-Keldysh effect. Coming from high fields, first the Wannier- 
Stark peaks are modulated by the Franz-Keldysh oscillations. However, 
as soon as the separation eFd between neighboring peaks becomes smaller 
than their spectral widths, the peaks can no longer be resolved individually 
so that only the Franz-Keldysh structure remains. 
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Figure captions 



Figure 1 Schematic illustration of the field-induced coherent motion of an elec- 
tronic wavepacket initially created at the bottom of a miniband. Here, the width of 
the miniband exceeds the LO-phonon energy Elo, so that LO-phonon scattering 
is possible. After Ref. UC 



Figure 2 Schematic representation of the transitions from the valence to the 
conduction band of a superlattice in the Wannier-Stark localization regime. After 
Ref g/. 

Figure 3 Full Bloch- oscillation dynamics corresponding to a laser photoexcita- 
tion resonant with the first-miniband exciton. (a) Time evolution of the electron 
distribution as a function of kp (b) Average kinetic energy, (c) current, and (d) 
THz-signal corresponding to the Bloch oscillations in (a). 

Figure 4 (a) Total THz-signals for eight different spectral positions of the ex- 
citing laser pulse: 1540, 1560, . . . , 1680 meV (from bottom to top), (b) Individual 
THz-signal of the electrons and holes in the different bands for a central spectral 
position of the laser pulse of 1640 meV. After Ref. 

Figure 5 (a) Total THz-radiation as a function of time; (b) Incoherently-summed 
polarization as a function of time. After Ref. 

Figure 6 (a) Electron and (b) hole contributions to the total THz-radiation of 
Fig. %a). After Ref 

Figure 7 Absorption spectra for various static applied electric fields for a 
GaAs/ Alo.sGao.j As superlattice (well (barrier) width 95 (lb) A). The vertical 
displacements between any two spectra is proportional to the difference of the 
corresponding fields. The Wannier-Stark transitions are labeled by numbers, the 
lower (higher) edge of the combined miniband by Ea (E\ ). After Ref. Ma/. 
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